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ABSTRACT 

In this paper, we construct full amplitude non- linear hydrodynamical models of fun- 
damental mode Galactic Cepheids and analyze the resulting theoretical period-colour 
and amplitude-colour relations at maximum, mean and minimum light. These theo- 
retical relations match the general form of the observed relations well. This agreement 
is, to some extent, independent of the mass-luminosity relations used, pulsation code, 
numerical techniques, details of the input physics and methods to convert theoretical 
quantities, such as bolometric luminosity and temperature, to observational quantities, 
such as V band magnitudes or (V — I) colours. We show that the period-colour and 
amplitude-colour properties of fundamental mode Galactic Cepheids with periods such 
that log(P) > 0.8 can be explained by a simple application of the Stefan-Boltzmann 
law and the interaction of the photosphere with the hydrogen ionization front. We 
discuss the implications of our results for explaining the behavior of Galactic Cepheid 
period-colour, and period-luminosity relations at mean light. 
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1 INTRODUCTION 



iKanbur fc Ngeowl i2004l paper I, hereafter KN) presented 
new observational characteristics of fundamental mode 
Cepheids obtained from an analysis of the Galactic Cepheid 
data and the OGLE LMC/SMC Cepheid data. These were 
the period-colour (PC) and the amplitude-colour (AC) rela- 
tions at maximum, mean and minimum light. KN analyzed 
these PC and AC diagrams in t he context of the work of 
ISimon. Kanbur fc Mihalasl Jl993t hereafter SKM), who used 
the Stefan-Boltzmann law and the fact that radius fluctua- 
tions are small in Cepheids, to derive the following equation, 
valid for optical pulsations of Cepheid variables: 



logT„ 



where T, 



10 



(14, 



V m 



(1) 



max I rain 



is the photospheric temperature at maxi- 



mum or minimum light. Consequently, if for some reason the 
PC relation is flat at maximum light, and the colour used is 
a good predictor of temperature, then equation (1) predicts 
a relationship between V band amplitude and T m i n , and 
thus a correlation between the V band amplitude and the 
observed colour (after extinction correction) at minimum V 
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band light. ICodel dl947l) found that Galactic Cepheids ex- 
hibit a spectral type that is independent of period at max- 
imum light. SKM analyzed existing Galactic Cepheid data 
to show that higher amplitude Galactic Cepheids are in- 
deed driven to cooler temperatures, and thus redder (B — V) 
colours, at minimum light in acc ordan c e with equation (1) 
and the observational findings of ICodel Jl94/f) . 

The reason why Galactic Cepheids follow a spectral 
type that is independent of period at maximum light was 
explained in SKM: it is due to the interaction of the pho- 
tosphere with the hydrogen ionization front (HIF). At max- 
imum light, the HIF is so far out in the mass distribution 
that the photosphere, taken to be optical depth 2/3, oc- 
curs right at the base of the HIF. Together with the HIF 
is a sharp rise in opacity. At this point, the mean free path 
goes to zero. In the absence of any significant driving, even 
though the surrounding atmosphere has a non-zero inward 
velocity, this opacity wall prevents the photosphere from go- 
ing deeper in the star and erases any "memory" of global 
stellar conditions. Thus for a large range of periods, the pho- 
tosphere occurs at the base of the HIF at maximum light at 
a temperature which is independent of period. At maximum 
light, this leads to a flat relation between period & temper- 
ature, period & (B — V) colour and period & spectral type, 
as seen in SKM. At other phases, the HIF lies too far in- 
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side the mass distribution to interact with the photosphere. 
This, together with equation (1), implies a relation between 
amplitude and (B — V) colour or temperature at minimum 
light. Equation (1) suggests that if T m i„ obeyed a flat re- 
lation with period, then there should be a relation between 
amplitude and colour at maximum light. Hence PC relations 
are connected to AC relations, and changes in the PC re- 
lations should be reflected in corresponding changes in the 
AC relations. 

Because of the extensive data now available, KN ana- 
lyzed recent Galactic and Magellanic Cloud Cepheid data 
in terms of PC and AC diagrams at maximum, mean and 
minimum light. Using (V — I) colours, they performed a 
more detailed analysis of the Galactic Cepheid data and 
found that the Galactic PC relation at maximum light dis- 
played a statistically significant break at 10 days, but was 
consistent with a single line at mean and minimum light. In 
the LMC, the PC relation displayed this break at all three 
phases, though the statistical significance at maximum light 
was marginal. The SMC PC relation displayed similar prop- 
erties to that in the Galaxy. 

Analysis of the Galactic Cepheid data in terms of AC 
diagrams confirmed the work of SKM, and extended it to 
the case of (V — I) colour: there is a relation between am- 
plitude and (V — I) colour at minimum V band light. For 
the LMC, short period Cepheids (log(P) < 1.0) show no 
relation between V band amplitude and (V — I) colour at 
minimum light, but long period Cepheids (log(P) > 1.0) are 
such that higher amplitude stars are driven to redder (V — I) 
colours and hence cooler temperatures at minimum light. At 
maximum light, short period Cepheids are such that higher 
amplitude stars are driven to bluer (V~ I) colours and hence 
hotter temperatures but long period Cepheids do not show 
such a relation. 

An understanding of the PC and AC relations derived 
by KN are important, not only for stellar pulsation and evo- 
lution studies of Cepheids, but also because the Cepheid 
period-luminos ity (PL) relation depends on the PC rela- 
tion (see, e.g., [Madore fc Freedmanlll99lh . The PL relation 
will reflect significant changes in the PC relation. Hence in 
studying PC and AC relations in different galaxies, we are 
studying the universality of the Cep h eid P L relation. Re- 
cent work bv iTammann et al.l (120021. l2003f). iFouaue et alJ 
j2003ft.lKanbur fc Ngeowl j2004l).lNgeow fc Kanburl i2004aTl ' 
ISandaee et alJ i2004l) and lstorm et alJ i2004| ) have sug- 
gested that the PL relation in the Galaxy is significantly 
different from that in the LMC and, further, that the 
PL relation in the LMC is non-linear. In contrast, cur- 
rent obser vations indicate that t he Galactic PL relation 
is linear iTammann et al] 120031 : iKanbur fc Neeowl 12004 
iNeeow fc Kanburll2004air ^ 

The purpose of this and subsequent papers is to con- 
front these new observed characteristics of Cepheids with 
the latest stellar pulsation models and interpret the results 
in terms of the theory presented in SKM, which has been 
summarized above. This approach will ultimately yield a 
qualitatively deeper understanding some of the reasons be- 
hind the variation of the Cepheid PL relation from galaxy to 
galaxy. Our ultimate goal is to use our theoretical models to 
estimate the effect of this variation quantitatively. For this 
paper we concentrate on Galactic Cepheid models which we 
compare with the same Galactic Cepheid data used in KN. 



The LMC/SMC Cepheid models will be presented in the 
forthcoming papers in this series. 

Our models and methods improve upon those in SKM 
in the following respects: First of all, they contain a for- 
mulation to model time dependent t urbulent convection 
jYecko et alJll99St IKoUath et al] 12002). in contrast to the 
purely radiative models used in SKM. In addition, we con- 
struct more models so that we can examine in greater de- 
tail the PC and AC characteristics of Cepheids. Secondly, 
we investigate the pulsation properties as a function of 
other phase points. When investigating PC and AC rela- 
tions at mean light, KN defined mean (V — I) colour as 
Vmean — Iphmean where V m ean is the V band magnitude 
closest to the value of Ao, the mean magnitude obtained 
from a Fourier decomposition, and Iphmean is the I band 
magnitude at this same phase. The reason for doing this is 
because the colour of the star at a certain phase is precisely 
the V band magnitude minus the I band magnitude at that 
same phase. In contrast, a definition of mean colour such as 
< V > — < I > folds in the phase difference (albeit small) 
that exists between the V and I band light curves. We adopt 
the definitions given in KN for colour at maximum, mean 
and minimum light. However, there are two phases when 
the V band magnitude is closest to the value of Ao: once 
on the ascending branch and once on the descending branch 
of the light curve. The colour (or the temperature), in both 
models and theory, need not necessary be the same at these 
two phases. In KN's analysis of the observed data, the de- 
scending branch mean was adopted as the "mean" . We pay 
specific attention to this detail in our results and discussion 
section. Finally, we look at (V — /) colours whereas SKM 
studied predominantly (B — V) col ours. The (V — I) colour 
is a good indicator of tempe rature iGonzalez fc Wallersteirj 
Il996t ITammann et alj feoOS) and is a crucial colour for the 
existing calibration of the extra-galactic distance scale (e.g., 
iFreedman et al]l200ll) . 



2 MODELS AND CODE DESCRIPTION 

The numerical techn iques and physics incl uded in these 
model s are detailed in lYecko et al] lll999f) and lKollath et all 
( 2002). The Florida code has been use d successfully t o mode l 



double mode Cepheid pulsationsjjKollaH^j^d] 



1998 



2000), 



Galactic first overtone Cepheids ]Feuchtinger_et_aL[ ' 
Ceph eid mass-luminosity (ML) relations <lBeluiheu^t"aLl 
2001) and the stu dy of Magellanic Clouds Cepheids 
feuchler et al. 1 120041) . The code takes the mass, luminosity, 
effective temperature, hydrogen and metallicity abundance 
(by mass), M, L, T e ,X, Z, as input parameters. In this study, 
we used (X,Z) = (0.70,0.02) as a representation of Galac- 
tic hydrogen and metallicity abundance. For the remain- 
ing parameters, mass and luminosity are connected by an 
adopted ML relation, hence it is only necessary to choose 
the mass and the effective temperature for computing a 
Cepheid model. In this study, we used two ML relations, 
which are calculated from evolutionary models appropriate 
for intermediate-mass stars: 



(i) The M L relation given b y IChiosil 1119891) . which was 
also used by ISimon fc Kanbuil dl995T) in their theoretical 
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Table 1. Input parameters for Galactic Cepheid models with pe- 
riods obtained from a linear analysis. The periods, Pq and Pi, 
are referred to the fundamental and first overtone periods, re- 
spectively. Similarly for the growth rate, rj. Both of the mass and 
luminosity are in Solar unit, the temperature is in K and the 
period is in days. 



M 


log(L) 




Po 




'/(> 


Pi 


m 




ML Relation from Bono et al. 


(2000) 




11.5 


4.279 


4830 


43.3924 


0. 


.042 


27.58 


-0.164 


11.1 


4.228 


4900 


37.8090 





.039 


24.40 


-0.157 


10.8 


4.188 


4950 


34.4156 





.037 


22.25 


-0.133 


10.5 


4.147 


5000 


30.8189 





.035 


20.26 


-0.119 


9.55 


4.009 


5075 


23.4835 





.025 


15.69 


-0.094 


9.45 


3.994 


5265 


20.0414 





.028 


13.55 


-0.051 


8.60 


3.856 


5300 


15.8313 





.022 


10.78 


-0.040 


7.70 


3.696 


5332 


12.1076 





.015 


8.291 


-0.036 


7.30 


3.618 


5440 


10.0124 





.014 


6.886 


-0.017 


7.30 


3.618 


5490 


9.69213 





.014 


6.670 


-0.007 


7.00 


3.557 


5490 


8.83394 





.013 


6.087 


-0.010 


6.80 


3.515 


5485 


8.31328 





Oil 


5.733 


-0.014 


6.45 


3.438 


5545 


7.12090 





Oil 


4.921 


-0.007 


6.10 


3.357 


5580 


6.16131 





.009 


4.267 


-0.006 


6.00 


3.333 


5590 


5.59661 





.006 


3.884 


-0.007 


ML Relation from Chiosi (1989) 


10.0 


4.534 


5150 


66.9094 


0. 


.144 


38.47 


-0.413 


8.50 


4.279 


5100 


44.8456 





.129 


27.38 


-0.238 


8.00 


4.184 


5090 


38.5673 





.116 


23.98 


-0.196 


7.40 


4.061 


5050 


32.4902 





.096 


20.55 


-0.171 


7.20 


4.019 


5250 


26.1109 





.081 


16.99 


-0.115 


7.00 


3.975 


5400 


21.8715 





.060 


14.45 


-0.077 


6.30 


3.810 


5390 


16.9148 





.055 


11.37 


-0.048 


5.50 


3.597 


5350 


12.4229 





.038 


8.445 


-0.051 


5.10 


3.478 


5396 


10.0049 





.030 


6.851 


-0.041 


5.00 


3.447 


5420 


9.37928 





.028 


6.436 


-0.036 


4.90 


3.416 


5440 


8.80955 





.026 


6.056 


-0.033 


4.80 


3.383 


5470 


8.21323 





.025 


5.657 


-0.027 


4.50 


3.282 


5460 


7.06755 





.017 


4.878 


-0.037 


4.40 


3.247 


5490 


6.56235 





.016 


4.536 


-0.032 


4.20 


3.174 


5560 


5.60322 





.015 


3.887 


-0.022 






Input parameters 


from SKM 




4.57 


3.306 


5707 


10.5472 


0.035 


7.234 


-0.000 


5.44 


3.578 


5550 


6.25949 


0.022 


4.340 


0.020 



Table 2. Temperatures at maximum and minimum light from 
full-amplitude non-linear model calculations. The periods, lumi- 
nosity and temperature are in days, Lq and K, respectively. 



Lmax Tmax L m l n T^min 



ML Relation from Bono et al. (2000) 



43.3924 


23323.69 


5301.44 


15265.45 


4653 


.25 


37.8090 


20547.96 


5355.92 


13844.78 


4759 


.77 


34.4156 


18675.64 


5399.02 


12811.00 


4681 


.60 


30.8189 


17046.03 


5446.49 


11592.02 


4719. 


.65 


23.4835 


12259.96 


5456.64 


8535.450 


4820 


.82 


20.0414 


12503.63 


5720.76 


7897.603 


4948 


.82 


15.8313 


8828.696 


5664.44 


6018.911 


5025 


.59 


12.1076 


5616.756 


5525.95 


4402.925 


5135. 


.11 


10.0124 


4470.551 


5531.94 


3703.182 


5274. 


.68 


9.69213 


4484.612 


5585.30 


3647.089 


5303 


.57 


8.83394 


3872.188 


5619.60 


3172.488 


5312 


.02 


8.31328 


3501.343 


5612.30 


2906.126 


5322 


.39 


7.12090 


2954.314 


5673.40 


2436.595 


5396 


.37 


6.16131 


2434.291 


5700.83 


2059.419 


5459. 


.73 


5.59661 


2260.885 


5691.47 


2023.329 


5513. 


.33 


ML Relation from Chiosi (1989) 


66.9094 


39933.45 


5468.07 


25451.85 


4807 


.32 


44.8456 


24179.71 


5738.10 


12478.80 


4707. 


.13 


38.5673 


19927.26 


5762.38 


9964.144 


4696 


.57 


32.4902 


15175.07 


5724.16 


7818.562 


4720 


.46 


26.1109 


13968.71 


5886.81 


7294.068 


4949 


.57 


21.8715 


13020.35 


6056.17 


7076.027 


5176. 


.93 


16.9148 


8608.982 


5999.40 


4782.056 


4961. 


.41 


12.4229 


5070.216 


5846.81 


3081.266 


4998 


.02 


10.0049 


3724.848 


5785.06 


2465.684 


5074. 


.73 


9.37928 


3427.493 


5770.41 


2320.390 


5115 


.94 


8.80955 


3130.614 


5748.80 


2187.659 


5156. 


.00 


8.21323 


2843.303 


5729.82 


2062.474 


5217. 


.61 


7.06755 


2080.732 


5572.38 


1722.053 


5295 


.72 


6.56235 


1885.999 


5565.48 


1587.056 


5332. 


.63 


5.60322 


1580.301 


5697.21 


1315.624 


5383. 


.56 


Input parameters from SKM 


10.5472 


2276.417 


5894.64 


1660.255 


5430.83 


6.25949 


5084.652 


6068.32 


2974.031 


5170.89 



study of long period Cepheids: 

log L = 3.61 log M + 0.924. (2) 

(ii) The ML relations given bv lBono et"aH (2000): 

logL = 0.90 + 3.35 log M + 1.36 logY -0.34 log Z, (3) 
= 3.35 log M + 0.726. 

The units for both M and L are in Solar unit. Note that 
equation (2) gives a much larger L/M ratio than equation 
(3). The range of L/M ratio covered by equations (2) and 
(3) is reasonably broad and encompasses a fair selection of 
ML relations in the literature. Our aim for this paper is not 
to test one ML relation against another, but rather, to show 
that the physical effect we are interested in is, to a large 
extent, independent of the ML relation chosen. 



The effective temperature, T e , was chosen to ensure a 
good fundamental growth rate and a stable fi rst overtone 
mode . After a linear non-adiabatic analysis iYecko et alJ 
1999) which yields the normal mode spectrum, an initial per- 
turbation, which is scaled to match the fundamental mode 
linear velocity eigenvector, is applied and the model followed 
until a stable limit cycle is reached. In some cases, this means 
following the pulsations for about 1000 cycles. This results 
in a full amplitude Cepheid oscillating in the fundamental 
mode. Following SKM, we can compute the temperature at 
an optical depth r = 2/3 and also the temperature and 
opacity profiles as a function of depth at various phases of 
the pulsation. 

In addition to the two ML relations used, we have also 
computed two models using the parameters (M, L, T e ) from 
table 1 of SKM, in order to compare the results based on 
purely radiative and the turbulent convective models, and 
further broaden the range of ML relation considered. The 
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input parameters for these models and for the two ML rela- 
tions used are presented in Tabled 

We compare theoretical quantities to observable quan- 
tities using a number of prescriptions: 

(i) We convert observed, extinction corr ected (V — I) to 
(B — V ) colours using the formula given in lTammann et al] 
( 2003), and then convert the (B — V) index to temperatures 
using the formula given in SKM. 

(ii) We use the BaSeL atmosphere database 1 llLeieund 
l2002l:IWestera et alJl2002l) to construct a fit giving tempera- 
ture as a function of (V—I) colour, with and without a log(g) 
term. Where a log(g) term is used, we obtain this, at the ap- 
propriate phase, from the models . This d atab ase is also used 
in, fo r example, iBeaulieu et al] (2001) and ICordier et al] 
(|2003j) for the study of Magellanic Cloud Cepheids. 

(iii) These techniques were also used to obtain the 
bolometric corrections required to convert theoretical lu- 
minosity curves into theoretical V band light curves. 
ISimon fc Kanburl lll995l) have shown that these corrections 
are small. 

The qualitative thrust of our results does not change ac- 
cording to the details of the method used to convert colours 
to temperatures and luminosity to V band light. In what 
follows, we report on results using method (ii) above with 
a log(g) term. One caveat here is that we use hydrostatic 
equilibrium atmospheres. The effective temperature and ef- 
fective gravity for these atmospheres are obtained from the 
relevant quantities at the photosphere in the models. For the 
effective gravity we use 

GM _ du 

where du/dt is again given by the instantaneous structure 
of the models. Such an approach has been widely used in 
the literature in recent years for both Cepheids and RR 
Lyraes (see, for examples, SKM, |Si mari_fc_i£anbjrJ I1995L 
Bono et al]ll999L ISandaee et aljll999l IBeaulieu et aljhoOll 
Caputo et all2002l . lKovacsi2003LlRuoDDO et all2004l) . Some 
justifi cation based o n d etailed modeling may be foun d in 
iKollath et all i2000T) and lKeller fc Mutschlecneil (Il970l) . We 
note that at most phases the du/dt term is small but can 
become large at certain phases during the ascending and 
descending branches. For the models considered here, this 
occurs for three models with log(P) between 1.5 and 1.65 at 
phases close to maximum light. 



3 RESULTS 

The results from the model calculations are summarized in 
Table Q |21 and |3] TableQgives the input parameters, as well 
as the computed periods and growth rates from the linear 
non-adiabatic analysis for the fundamental and first over- 
tone modes for the models constructed in this paper. We 
note that the models used generally had reasonably high 
fundamental mode growth rates and negative first overtone 



1 http : //www. astro .mat .uc .pt/BaSeL/ 



growth rates. This puts our models well inside the funda- 
mental mode instability strip for all periods considered. Ta- 
ble |21 gives the temperatures at the maximum/minimum 
luminosity at full amplitude for the corresponding models 
in Table Table |3] gives similar quantities but for the as- 
cending and descending branch means. In Table |3] < L > 
is the average luminosity obtained from the light curve in 
one pulsating cycle. L mea „ is the luminosity that is ob- 
tained from the closest temporal grid points (from the light 
curve) to the < L >, on both of the ascending and descend- 
ing branches. T mean are the corresponding temperatures at 
these grid points. Since < L > and L mean may not be equal 
to each other, we obtain the temperatures that correspond 
to < L > by interpolation between L mea „ and the adjacent 
luminosity that bracket the < L >. These temperatures are 
referred as T^^n in Table |3] Note that the temperatures 
at mean light on the descending branch is normally cooler 
than the mean light temperatures on the ascending branch. 

3.1 The PC relations 

In this subsection we compare the PC relations from the 
model calculations and the observations. Figure shows a 
four panel plot of log-period against extinction corrected 
(V — /) colour for Galactic Cepheids as given in KN. Su- 
perimposed on this are values obtained from our models, 
where the theoretical temperatures are converted to (V — /) 
colours using method (ii) described above. The sol i d and 
ope n squares represent models based on the IChiosl ]l989T) 
and iBono et al] <|2000T) ML relations respectively, whereas 
the solid circles refer to models computed with the ML re- 
lation used in SKM. The left panels are for maximum (top) 
and minimum (bottom) light respectively, whilst the right 
panels are for the means of the ascending and descending 
branches, respectively. We see that there is very good agree- 
ment between models and observations in this diagram. Fig- 
ure |21 shows the same quantities but on the log(P)-log(T) 
plane. The observed (V — /) colours for the Galactic Cepheid 
data are converted to the temperature by using the BaSeL 
database mentioned in previous section, where the log(<?) 
terms for the observed colo urs are approximated by using 
log(p) = 2.62 - 1.211og(P) iBeaulieu et aljl200lD . We note 
that the qualitative nature of our results is the same in Fig- 
ure and |21 the flat nature of the relation at maximum 
light, the non-zero slope at minimum light and an amalgam 
of these at mean light. 

Recall that KN analyzed a large sample of Galactic 
Cepheid data and showed that the PC relation was con- 
sistent with a single line at mean and minimum light, but at 
maximum light there was significant statistical evidence for 
a break at 10 days. Figures & [2] suggests that this break 
actually occurs at log(P) ~ 0.8. As in KN, we can perform 
a two line regression, for Cepheids with periods smaller and 
greater than log(P) = 0.8, and compare with a one line re- 
gression across the entire period range. This is what is meant 
by short and long period in the following paragraph. Table^l 
shows the phase, overall dispersion, slope and error at that 
phase and then these same quantities for the short and long 
period Cepheids respectively. From Figure |2| and Table |11 
we can conclude the following results: 

(i) The overall slope in the log(P)-log(T) plane at maxi- 
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Figu re 1. The period-c olour (PC) relations for the Galactic Cepheids. The crosses are the observed d ata points, which are taken from 
iKanbur fc Ngeowl 12004) . The open and solid squares are the Galactic models with lBono et all 1200(11) and IChiosil Il989l) ML relations, 
respectively, whilst the solid circles are results for models computed with the ML relation used in SKM. The temperatures from the 
models are converted to the (V — I) colour using the BaSeL database with the log(rj) term. Left (a): PC relations at maximum and 
minimum light; Right (b): PC relation at mean light, for both of the acceding and descending branches. 
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Figure 2. The plots of log(T) vs. log(P) for the Galactic data and models. The symbols are same as in Figure Q] The conversion from 
observed (V — I) colour to temperature is done with the BaSeL database. Left (a): The log(P)-log(T) plots at maximum and minimum 
light; Right (b): The log(P)-log(T) plots at mean light, for both of the ascending and descending branches. 



mum light is flat and is significantly different from the slope 
at mean and minimum light. 

(ii) The short period slopes in this same plane are signif- 
icantly different from zero at all three phases. In contrast, 
the long period slope at maximum light in this plane is close 



to zero and is significantly different from the slope at mean 
and minimum light. 

(iii) If we reject all Cepheids with log(P) < 0.8 and test, 
using the methodology of KN, whether the remaining data 
are consistent with a single line or with two lines broken 
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Figure 3. The V band amplitude-colour (AC) relations for the Galactic Cepheids. The symbols are same as in Figure^ The bolometric 
model light curves are converted to V band light curves with the bolometric corrections obtained from BaSeL database. Left (a): AC 
relations at maximum and minimum light; Right (b): AC relation at mean light, for both of the acceding and descending branches. . 



at a period of 10 days, we find that in the log(P)-rog(T) 
plane, even at maximum light, the data are consistent with 
a single line. Thus the broken PC relation at maximum light 
for Galactic Cepheids described by KN is due primarily to 
the different PC relation obeyed by short (log(P) < 0.8) 
period, fundamental mode Cepheids. 

We note that this good agreement exists for maxi- 
mum, mean and minimum light. The exceptio n is perhaps 
the lo nger period models constructed with the lBono et alJ 
(2000) ML relation, with log(P) > 1.3, at maximum light. 
We discuss these models in Section 4. Since the observed 
data exhibit very similar characteristics in both period- 
color and period-temperature planes, that is, a distinctly 
non-linear relation with a sudden break at a period of 
0.75 < log(P) < 1.0 and since the models agree well with 
the data, it follows that existing models contain the physics 
to explain these observations. We study this in Section 3.3 

3.2 The AC relations 

Figure [3] depicts plots of the AC relations for both models 
and data, at maximum, minimum, ascending and descend- 
ing branch means. The bolometric corrections obtained from 
the BaSeL database are used to construct the V band ampli- 
tudes from the model light curves. We see good agreement 
between models and data in these di agrams, though m odels 
constructed with our code using the iBono et alJ (l200Cl) ML 
relation have smaller amplitudes even though these mod- 
els match observa tions in the AC diagrams. Models con- 
structed using the Chiosi ( 1989J) ML relation with our code 
generally have larger amplitudes. We caution here that the 
computed amplitudes, unlike the periods, depend strongly 
on the strength of the artificial viscosity and of the assumed 
turbulent eddy viscosity. 



Table gives the results from a linear regression of V 
band amplitude against log T at the three phases of interest 
for the Galactic Cepheid data, taken from Figure |5] & El 
Equation (1) shows that, if the variation of \ogT max with 
period has a shallow slope close to zero, then an equation of 
the form logT m i„ = a + bV a m P is such that b is about —0.10. 
From Table|KJ we see that the overall slope is —0.079, which 
is very close to the theoretically expected value. Accounting 
for the error on the slope shows that the actual slope is 
statistically indistinguishable from the theoretical value of 
—0.10. In addition, Table |S] implies that the overall slope 
for Vamp-log T m ax relation is close to zero across the entire 
period range. This strongly supports the validity of equation 
(1) and our interpretation of the predictions arising from this 
equation if either log T max or log T m in obey a flat relation 
with period. Figure 21 argues this graphically with a two 
panel plot displaying V band amplitude against \ogT m in 
and logTmax on the left and right panels, respectively. An 
AC relation may also occur if either log T max or log T m i n has 
a shallower relation with period than the other quantity. 



3.3 The interaction of the HIF and photosphere 

In Figures 15171 we present the plots of the temperature and 
opacity profiles, i.e. temperature and log opacity against lo- 
cation in the mass distribution respectively, for a representa- 
tive selection of the models in Tabled We display plots for a 
long period (log(P) > 1.0), short period (log(P) < 1.0) and 
a 10 day period model. The mass distribution is measured 
by the quantity Q — log(l — M r /M), where M r is the mass 
within radius r and M is the total mass. Each panel shows 
the temperature and opacity profiles at maximum, mean and 
minimum light with dotted, solid and dashed curves, respec- 
tively. For each model, we also include the plots for the case 
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Table 3. Temperatures at mean light from full-amplitude non-linear model calculations. See text for the meanings of < L >, L 
Tmean and Tmean- The periods, luminosity and temperature are in days, Lq and K, respectively. 
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4147.76 
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5563.73 
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2738.62 
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5623.67 


2739.4820 
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2271.74 
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5533.36 


5660.78 
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4867.35 


32.4902 


11522.9 


11331.742 


5312.26 


11511.286 


4796.29 


5335.61 


4797.57 


26.1109 


10436.0 


10650.701 


5571.95 


10432.457 


4983.82 


5541.43 


4984.29 
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5688.05 
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3007.79 
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5557.55 


5239.67 
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2417.35 
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5369.50 
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Input parameters from SKM 


10.5472 


2022.70 


2028.8310 


5837.90 


2031.3158 


5634.15 


5832.38 


5627.78 


6.25949 


3777.86 


3731.9089 


5710.18 
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5373.31 
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Table 4. log(T) vs log(P) regression at maximum, mean and minimum light for the Galactic Cepheid data, a is the dispersion from the 
regression of the form: log(T) = a + blog(P). The subscripts £ and g refer to long (log(P) > 0.8) and short period Cepheids, respectively. 



phase 


&all 


ball 


a all 




b L 




&S 


bs 


as 


Max 

Mean 

Min 


0.016 
0.017 
0.022 


-0.033 ±0.005 
-0.062 ± 0.006 
-0.077 ±0.007 


3.807 ± 0.005 
3.801 ± 0.005 
3.786 ± 0.007 


0.017 
0.019 
0.023 


-0.014 ±0.009 
-0.068 ±0.010 
-0.076 ±0.012 


3.787 ± 0.009 
3.808 ±0.010 
3.786 ±0.012 


0.014 
0.014 
0.019 


-0.109 ±0.023 
-0.131 ±0.023 
-0.193 ±0.031 


3.862 ± 0.016 
3.846 ± 0.016 
3.865 ± 0.022 



Table 5. V band amplitude (V am p) vs. logT m i„ regression for the Galactic Cepheid data, cr is the dispersion from the regression of the 
form: logT m i n = a + bVamp- The subscripts £, and g refer to long (log(P) > 0.8) and short period Cepheids, respectively. 



phase 


Call 


ball 


Call 




b L 




&S 


bs 


as 


Max 

Mean 

Min 


0.018 
0.021 
0.024 


0.011 ±0.007 
-0.050 ± 0.009 
-0.079 ±0.010 


3.770 ± 0.006 
3.784 ± 0.007 
3.778 ± 0.008 


0.017 
0.021 
0.023 


0.018 ±0.008 
-0.048 ± 0.009 
-0.073 ±0.011 


3.757 ± 0.007 
3.776 ± 0.009 
3.766 ± 0.009 


0.015 
0.017 
0.023 


0.038 ±0.015 
-0.007 ±0.017 
-0.040 ± 0.023 


3.760 ±0.011 

3.761 ± 0.012 
3.760 ±0.016 
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when the mean is on the ascending and descending branch 
respectively. The photosphere is marked with a filled circle 
in these figures. Figure|H]illustrates the luminosity and tem- 
perature profiles against mass distribution in one pulsating 
cycle. 

In general, the HIF sweeps back and forth in the mass 
distribution as the star pulsates. A naive way of thinking 
about this would be that the HIF is furthest out and furthest 
in, in the mass distribution, at maximum and minimum light 
respectively (SKM). Our models indicate that the situation 
is more complicated than this, and we summarize the main 
results from Figure 15171 as follow: 

(i) At minimum light, the photosphere is closest to the 
HIF at short period but moves progressively away from the 
base of the HIF as the model period increases. This is in 
concordance with the lowest panel of Figures Q and |2] which 
display the period temperature relation for models and data 
as one single line such that longer period models have cooler 
photospheric temperatures at minimum light. However, at 
maximum light, the photosphere lies at the base of the of 
the HIF for all the models presented. 

(ii) We can quantify this by defining the quantity A to be 
the "distance", in units of Q, between the photosphere and 
the point in the HIF where dT/dQ is a maximum. Figure [5] 
illustrates this and Figure [TU1 portrays the results when we 
plot A against log(P) at maximum, minimum and the as- 
cending and descending branch means. In Figure the er- 
ror bars are estimated from the coarseness of the grid points 
around the location of HIF. We see clearly that at maxi- 
mum light, A is constant over a large period range, whereas 
at minimum light, A gets larger as the period increases. 
The value of A at the phases of ascending and descending 
branch means is a mixture of the behavior at maximum and 
minimum light. There is an indication that for long period 
models (log(P) > 1.3) at maximum light, the photosphere 
starts to become disengaged from the base of the HIF. As the 
period increases, so does the L/M ratio and, in general, the 
effective temperatures become lower. Kanbur (1995) showed 
that such changes push the HIF further inside the mass dis- 
tribution, making it harder to interact with the photosphere 
even at maximum light when it is closer to the surface. This 
can cause a non-zero slope in the log(P)-log(T ma:c ) plane at 
long periods. This was also noted by SKM. 

(iii) The HIF at mean light is not always between its lo- 
cations at maximum or mean light. Figure |H| illustrates our 
point and shows various locations on the ascending and de- 
scending branch of a typical light curve from our models. 
The lower two panels show the temperature profile at each 
of these phases. We see that the HIF at phase b is actually 
further out in the mass distribution than at maximum light 
and conversely the HIF is further in in the mass distribution 
at the phase corresponding to mean light on the descending 
branch than at minimum light. Moreover, the nature and 
shape of the HIF seems to change significantly during the 
pulsation period. 

(iv) At short/long period, the HIF at ascend- 
ing/descending branch mean light is located in front /behind 
its location at maximum/minimum light. 

(v) We also note that the nature of the HIF changes 
when it is located outside of its locations at either maxi- 
mum or minimum light. These changes in the HIF occur 
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Figure 9. Illustration of the location of HIF and photosphere, 
and the way to calculate A. The crosses represent the zones in 
the model calculations. The location of the photosphere, at r = 
2/3, and the HIF are marked with filled circles. The location of 
the HIF is defined as the zone with the steepest gradient in the 
temperature profile. Therefore, A is the "distance" (in terms of 
log[l — M r /M\) between the photosphere and HIF. Upper and 
lower panels illustrate the cases that the photosphere is far away 
and close to the HIF, respectively. 

at period ranges corresponding to the location of bumps on 
the ascending or descending branches of the light curve: the 
Hertzsprung progression. 

(vi) If we define "amplitude of HIF oscillation" as the 
range of log(l — M r /M) values occupied by the HIF, then 
there is a correlation between this and the resulting ampli- 
tude of the photospheric light curve. 

3.4 Cepheids with log(P) < 0.8 

Figure and the results of Tables 01 & El suggest that 
Cepheids with log(P) < 0.8 do not obey a flat relation with 
period at maximum light. Our current grid of models does 
not include any models with such periods, but Figures 1101 
does indicate that the photosphere will be located at the 
base of the HIF at all phases of pulsation as is the c ase with 
RRab stars jKanbuJll995l : iKanbur fc Phillip jri996tl . Figure 
1101 also indicates that as the period gets shorter the "dis- 
tance" in the mass distribution between the HIF and the 
photosphere is the same at maximum and minimum light. 
This does not necessarily mean that there w ill be a flat PC 
rela tion at both these phases a s explained bv lKanbuil dl995ll 
and iKanbur fc Phillips! (1996) for the case of RRab stars. 
These short period, fundamental mode Cepheids will be in- 
vestigated in a future paper. 



4 LIGHT CURVE STRUCTURE 

It is desirable to compare the theoretical light curves from 
models quantitatively to the obser ved light curves, for ex- 
ample by Fourier decomposition JPav ne-Gap oshki r3 ll947l : 
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Figure 4. The plots of log(T) vs. V amplitudes (V am p) for the Galactic data and models. The symbols are same as in Figure The 
conversion from observed (V — I) colour to temperature is done with the BaSeL database (see text for details). Left (a): The plot of 
Vamp-log(T m i n ), The dashed line is not the best-fit regression to the data, but illustrates how closely the overall trend of the data 
follow a slope of —0.1, as predicted from equation (1). This line is constructed with log(T m i„) = 3.78 — O.lVamp- Right (b): The plot of 
Vam P -log(Tmax). The dashed line show the average value of log(T ma x), where < log(T ma x) >= 3.78. 




Figure 5. The temperature and opacity profiles for a long period model. The dotted, solid and dashed curves are for the profiles at 
maximum, mean and minimum light, respectively. The filled circles are the location of the photosphere at r = 2/3 for each phases. The 
mean light profiles at the ascending and descending branch are the solid curves that lie close to the profiles at maximum light (dotted 
curves) and minimum light (dashed curves), respectively. Left (a): Temperature profile; Right (b): Opacity profile. 
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Figure 6. Same as Figure l5l but for model with period at 10 days. 
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Figure 7. Same as Figure 151 but for model with short period. 



ISimon fc Kanbuilll995l) or Principal Component Analysis 
{Kanbur et alJl2002h . The three panels of Figure [TT] present 
the plots of the Fourier amplitudes Ai & A2 (top and bottom 
left-side panels respectively) and <j>2i — (j>2 — 20i (right-side 
panel) against log(P) for the observed data and the model 
light curves. The Galactic Cepheid data, which are the same 
as those used in KN, and the V band model light curves are 
subjected to a Fourier decomposition of the form: 



fc=JV 

V = Ap + ^ [A k cos{kut + <p k )], 

k-l 



(4) 



where ui = 2n/P, and P is the period. We used a simu- 
lated annealing technique to perform this Fourier decompo- 
sition. This method significantly reduces numerical noise in 
the Fourier decomp osition of sparsely sampled periodic data 
jNeeow et alJl2003ft . 

In terms of these diagrams, there is much better agree- 
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Figure 8. Illustrations of the temperature profiles at different pulsation phases, where q = M r /M . The upper and lower panels are the 
light curves and the temperature profiles at various phases, respectively. The lower-left panel shows the changes of temperature profiles 
from minimum light to maximum light (ascending branch), and the lower-right panel shows the similar temperature profiles in descending 
branch. In lower panels, the dashed curves represent the temperature profiles at maximum and minimum light, while the solid curves 
represent the temperature profiles at other phases as indicated in upper panels. 



merit between models and observations when we use the 
Chiosi ML relation (equation [2]). These diagrams also il- 
lustrate that the agreement between models and data in the 
log(P)-temperature or log(P)-colour plane is a necessary, 
but not sufficient, condition for agreement between observed 
and theoretical light curves. This also emphasizes that the 
physical effect of interest in this paper - the interaction of 
the photosphere and HIF - is a fundamental physical effect 
which is, to some extent, independent of the ML relation 
used. 

The Florida code used in this study contains a recipe 
to calculate the time dependent turbulent convection whilst 
SKM found a similar result using a purely radiative, dy- 
namically zoned code. Our results imply that this interac- 
tion is also independent of numerical techniques and physics 
included in the code. What is needed is the presence of hy- 
drogen and an accurate description of the opacity when hy- 
drogen starts to ionize. 

At first sight, the left panels of Figurc HTI might lead one 
to question the ML relation in equation (3). However, the 
numerical recipe to model convection co ntains seven dimen - 
sionless parameters (the a parameters, lYecko et al.lll999l) . 
many of which can influence the amplitude. A more de- 
tailed study and comparison of the effect of ML relations 
may remove this discrepancy. For example, the two models 
constructed with the ML relation used in SKM, shown as 
black dots in Figure fTTI do lie closer to the observations in 
the Fourier plane. This may also explain the discrepancy 
between observations and theor y in the log(P)-(V — I) m ax 
plane for log(P) > 1.3 when the lBono et alJ feOOfl ML re- 
lation is used. 



5 CONCLUSION AND DISCUSSION 

By looking at the way the Cepheid photosphere, the region 
where the Cepheid continuum is generated, interacts with 
the HIF, we have provided a simple qualitative physical ex- 
planation for the observed PC properties of fundamental 
mode Galactic Cepheids with log(P) > 0.8. This explana- 
tion relies on the fact that the opacity in a Cepheid becomes 
very high when hydrogen starts to ionize. This acts as a wall 
and prevents the photosphere going any deeper, and leads to 
a flat relation between period and temperature at the phase 
when the HIF interacts with the photosphere. For Galac- 
tic Cepheids this is observed at maximum light. At other 
phases, the photosphere is located away from the opacity 
wall and its temperature is related to the global properties 
of the star and hence its period. This explains, convincingly, 
why the Galactic Cepheid period-temperature relation is flat 
at maximum light and has a non-zero, single slope at mean 
and minimum light. 

The qualitative nature of this idea remains true whether 
the pulsation code is purely radiative (as in SKM), or has 
a numerical recipe to model time dependent turbulent con- 
vection as in this work. The interaction between the photo- 
sphere and HIF ma y also provide s o me ex planation for the 
suggestion made bv iKervella et all (I2004T) . that the region 
where spectral lines are formed do not necessary move ho- 
mologously with the region where the K-band continuum is 
formed. In addition, because we have used two ML relations 
with a wide range of L/M ratio in our study, and because the 
same physical effect is present in models constructed with 
either relation, the qualitative nature of our result is inde- 
pendent of the ML relation used. However, the ML relation 
used and its variation with metallicity, will dictate how the 
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Figure 10. The plots of A as function of log(P). The symbols are same as in Figure^ The dashed lines represent (r 
boundary of the HIF. 
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Figure 11. The plots of first two Fourier amplitudes, as defined in equation (4), and <p2l as function of period for the observational data 
and models. The symbols are same as in Figure HI Upper Left (a): Plot of A± vs log(P); Lower Left (b): Plot of A 2 vs log(P). Right (c): 
Plot of (f>21 vs log(P). 



interaction of the HIF with the photosphere changes with 
Cepheids in a different metallicity environment. In the next 
paper in this series, we will investigate how these ideas can 
be used to explain the PC, AC and PL properties of funda- 
mental mode LMC Cepheids. 

We have also found evidence that the non-linear na- 
ture of the Galactic PC relation at maximum light reported 
by KN is due to short period fundamental mode Cepheids 



with log(P) < 0.8. These short period stars follow distinctly 
different PC relations and deserve further detailed study. 



5.1 The effect on the PL relation 

What bearing do the results of this paper have on the 
Cepheid PC and PL relations at mean light? It is clear that 
PC relations at different phases contribute to the PC rela- 
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tion at mean light. If we choose (V — I)i as the extinction 
corrected colour at phase i, then a PC relation at this phase 
is: 

(y-I)i = ai + bi log(P), 

and taking the average over a pulsation period, i = 1...N, 
this becomes: 

i-N i=N 

<v >-</>= i$>* + ^X> log ( p )- ( 5 ) 

1=1 1=1 

It is clear that the average intercept and slope will be 
affected by their values at ind ividual phase points such 
as maximum or minimum light. iNeeow fc Kanbul J2004bh 
compute PC relations for Galactic and Magellanic Cloud 
Cepheids at all phases between and 1 and show that, 
for example for Galactic Cepheids, as the phase approaches 
maximum light, the slope of the PC relation becomes flatter. 
Since the mean ligh t PC relation a ffects the mean light PL 
relation jMadore fc Freedmarll99llL our aim of understand- 
ing the reasons behind changes in the PC and AC relations 
at different phases has a direct bearing on understanding at 
least one cause of the possible variation of the mean light 
Cepheid PL relation from galaxy to galaxy. Below we out- 
line how our results are pertinent to studies of the variability 
of the Cepheid PL relation from galaxy to galaxy, though 
we emphasize that much of this discussion is dependent on 
theoretical model and data analysis currently being under- 
taken. Hence the following discussion is a road-map of some 
aspects of our future work. 

The evidence that shows the slope of the LMC 
PL relation at mean light is significantly different 
to the slope of the mean lig ht PL relation in th e 
Galaxy has been prov ided in iTam mann et alJ i2003l) . 
INeeow fc Kanburl J2004al) and lStorm et alJ J2004h Further- 
more, the mean light LMC PL relation suffers a change 
at a period of 10 days whilst the mean light Galac- 
tic PL relation dat a is consistent with a single slope 
with current data |Ta mmann et alj 120031: IFqu que et all 
| 2003l: IKanbur fc Ngeowl 12004 i Ngeow fc Kanburl l2004af 
bandage et al.ll2004l ; Istorm et al.lho04l) . 

To understand the effect of the PC relation on the PL 
relation, consider the period-luminosity-color (PLC) relation 
valid at any phase: 

M v =a + 61og(P)+c(V-7). 

If the PC relation is broken, say at a period Po, then we 
have: 

(V -I)=x + ylog(P),P <P . 
(V -I)=x' + y' log(P) I P>P . 

Substituting these two equations into the PLC relation leads 
to two PL relations: 

M v = a + cx + (b + cy) log(P), P < P Q . 

M v = a + cx' + (b + cy') log(P), P > P . 

Hence changes in the slope and intercept of the PC relation 
have a direct bearing on the slope and intercept of the PL 
relation . However, the results from this paper and also from 
KN and lNeeow fc Kanburl (l2004bD imply that the mean light 



PC relation is affected by changes in the PC relation at dif- 
ferent phases during a pulsation cycle, as given in equation 
(5). Thus the study of PC relations at various phases im- 
pacts on the variability of the mean light PL relation from 
galaxy to galaxy. 

In this paper, we have updated and extended the work 
of SKM to provide an account of a simple physical mecha- 
nism, the interaction of the photosphere and HIF, which can 
change the properties of PC relations for Galactic Cepheids 
and so affect the PL relation. Specifically, a new result of 
our work is that for Galactic Cepheids with log(P) > 0.8, 
the maximum light PC relation is flat, i.e., the HIF and 
photosphere are engaged at maximum light. We also study 
the changes in AC relations because they are linked to PC 
relations through equation (1). A flat PC relation at maxi- 
mum light leads to an AC relation at minimum light of slope 
~ — 0.1 from equation (1), as seen in Figure [I] and Table [K] 

However, a Cepheid in the LMC will have a different ML 
relation, usually in the sense that LMC ML relations have 
higher L/M ratios. Then in order to compare Galactic and 
LMC Cepheids of t he same period, the L MC Cepheid needs 
to be hotter CKN). iKanburl (119951) and IKanbur fc Phillip^ 
lll996ft found that such changes in the (M,L,T e ) triplet 
describing the model changes the relative location of the 
HIF and hence the phase and the range of periods at which 
they interact. Our hypothesis is that for LMC models, it 
is only after log(P) > 1.0 that the HIF and photosphere 
are engaged at phases around maximum light, leading to 
a flat PC relation. This flat PC relation for long period 
(log(P) > 1.0) LMC Cepheids is one cause for the non- 
linear na ture of the mean light LM C PL relation (KN). Em- 
pirically, INeeow fc Kanburl i2004bh provide preliminary evi- 
dence that the LMC PC relation is flat only for log(P) > 1.0 
whereas this crossover period is log(P) ~ 0.8 in the case of 
the Galaxy. In case of SMC Cepheids, KN also provide evi- 
dence that the maximum light PC relation is not flat even for 
long period Cepheids. Our contention is, using some of the 
results of this paper, that amplitudes in SMC Cepheids are 
not high enough to force an interaction between the photo- 
sphere and HIF , as the SMC Cepheids have lower amplitudes 
than the LMC JPaczvriskv fc PindoJ200fj) . This again must 
wait confirmation from a new set of SMC models. 

The final effect on PL relations at mean light will de- 
pend on the behavior at other phases, such as minimum 
light. For example, preliminary calculations to be presented 
in a future paper imply that certain changes seen in the PC 
relation at minimum light (KN) may correspond to the first 
overtone mode becoming stable. Not surprisingly, a thor- 
ough quantitative study of this must await the analysis of 
LMC and SMC Cepheid models, which are currently under 
construction. The details of the influence of flat/non-fiat PC 
relations on the mean PC and PL relations, the crossover pe- 
riod of Po that could vary from galaxy to galaxy, and the 
comparison of PL relations at maximum and minimum light 
for different galaxies are beyond the scope of this paper, and 
will be presented in the future papers. 

iGroenewegen et alJ i2004T) recently reported a metallic- 
ity dependence in the zero point of Cepheid PL relation at 
mean light. They used 34 Galactic Cepheids with individual 
metallicity measurements and then supplemented this sam- 
ple with primarily long period Magellanic Cloud Cepheids 
to show the existence of a quadratic term in log(P) in the 
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PL relation. When they used primarily Galactic Cepheids 
in their sample, they found no evidence of a quadratic term. 
They interpreted these results as being due to a metallic- 
ity dependent zero point in the PL relation. These results 
are also consistent with the Cepheid LMC PL relation hav- 
ing different slopes fo r long and short per io d Cepheids as 
sugges ted by KN and lSandaee et all i2004) . ISandage et alJ 
(2004) also plotted amplitude-mean colour relations in lu- 
minosity bins. However, KN's AC relations were along the 
instability strip. It will be interesting to apply the precepts 
behind equation (1) and multi-phase AC relations in lumi- 
nosity and/or period bins, that is, across the instability strip. 
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